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地 磁场 磁力 线 可 视 化 种 子 点 选取 的 
磁场 强度 线 积分 等 分 算法 
高 广大 大 晓 序 ”名 佳 ! 分 自明 : 


1( 中 国 科学 院 国家 空间 科学 中 心 ” 北 京 100190) 
2( 中 国 科学 院 大 学 ”北京 100049) 


摘 要 将 磁力 线 按 流 线 进行 绘制 是 地 磁场 可 视 化 的 通用 方法 ， 磁 力 线 种 子 点 的 选取 决定 了 绘制 磁力 线 的 玻 密 ， 
其 蔓 密 程度 能 和 否 真实 反映 地 磁场 强度 分 布 是 评价 地 磁场 可 视 化 效果 的 关键 . 基于 磁 经 圈 均 匀 角 度 种 子 点 选取 算法 
绘制 的 磁力 线 通常 不 能 客观 反映 地 磁场 强度 的 空间 分 布 , 针对 这 一 不 足 , 提出 一 种 等 分 磁场 强度 线 积分 的 磁力 线 种 
子 点 选取 算法 . 利用 该 算法 对 地 磁场 IGRF 模型 和 工 96 模型 描述 的 地 磁 内 外 源 场 进行 可 视 化 绘制 ,对 磁力 线 追 踪 
结果 中 出 现 的 元 余 磁 力 线 进 行 过 滤 , 统计 分 析 了 绘制 磁力 线 的 空间 分 布 与 地 磁场 强度 空间 分 布 的 相关 性 , 结果 表明 
该 算法 能 够 较 好 地 实现 对 地 磁场 的 可 视 化 . 

关键 词 地 磁场 可 视 化 , 种 子 点 选取 , 宛 余 磁力 线 过 滤 , 线 积分 等 分 算法 
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Abstract Drawing the magnetic lines as streamlines is a general method of visualizing geomagnetic 
field. A key factor to evaluate the effect of the geomagnetic field visualization is whether the space 
distribution of geomagnetic field lines is consistent with that of magnetic field intensity, while the dis- 
tribution of geomagnetic field lines is determined by seed point selection. The traditional algorithms 
that select seed points with uniform angles on magnetic meridian circles cannot objectively reflect 


the space distributions of magnetic field intensity. This paper proposes an algorithm of selecting 


* 中 国 科学 院 空间 科学 科技 领域 云 建 设 与 应 用 项 目 资助 (XXH12503-05-08) 
2015-03-18 收 到 原稿 , 2015-05-13 收 到 修 定稿 
E-mail: daoqlxQ@163.com. 通信 作者 邹 自 明 , E-mail: mzou@nssc.ac.cn 


148 


Chin. J. Space Sci.。 空间 科学 学 报 ”2016, 36(2) 


seed points with equal line integral of magnetic field intensity. The algorithm is applied to draw the 


geomagnetic field lines with the data from T96 model and IGRF model. The redundant magnetic 


field lines existing in the result are removed. Statistical analysis and comparison between the space 


distribution of magnetic field lines and the geomagnetic field intensity reveals that this algorithm 


can effectively visualize the geomagnetic field. 
Key words 


removing, Equal line integral algorithm 


0 引言 


地 磁场 可 视 化 可 以 直观 展现 地 磁场 的 形态 结构 
与 强度 分 布 等 物理 特征 , 通常 认为 地 磁场 是 一 种 典型 
的 矢量 场 , 将 地 磁场 磁力 线 按 流 线 进行 绘制 , 可 以 形 
象 描绘 磁场 结构 , 其 是 一 种 对 地 磁场 进行 可 视 化 的 通 
用 方法 . 

有 效 的 地 磁场 可 视 化 算法 应 能 够 表现 出 地 磁场 
的 物理 特征 . 在 磁场 强度 越 大 的 区 域 , 单位 面积 内 的 
磁 通 量 越 大 , 表现 在 视觉 效果 上 的 磁力 线 分 布 应 该 越 
密 ; 反之 , 在 磁场 强度 越 小 的 区 域 ， 则 磁力 线 应 该 越 
芯 . 因此 , 评价 地 磁场 磁力 线 可 视 化 效果 的 关键 因素 
是 绘制 磁力 线 的 朴 密 与 磁场 强度 空间 分 布 之 间 的 相 
关 性 . 利用 磁力 线 对 地 磁场 进行 可 视 化 时 , 对 选取 的 
每 个 磁力 线 种 子 点 进行 追踪 都 将 生成 一 条 磁力 线 , 所 
以 种 子 点 的 空间 分 布 直接 决定 了 磁力 线 的 空间 分 布 . 
选取 出 能 够 真实 反映 地 磁场 强度 空间 分 布 的 种 子 点 ， 
将 提高 地 磁场 磁力 线 可 视 化 的 有 效 性 . 

目前 对 地 磁场 磁力 线 可 视 化 的 工作 一 般 是 在 地 
球 磁 经 圈 、 磁 纬 圈 、 磁 赤道 面 或 菜 个 经 纬 网 格 区 域内 
选取 种 子 点 . 但 通常 的 做 法 是 选取 空间 均匀 分 布 ( 例 
如 均匀 角度 ) 的 种 子 点 , 缺乏 对 种 子 点 空间 分 布 的 物 
理 约束 . 例如 , Tsyganenkol 在 与 太阳 风 方 向 平行 的 
磁 经 圈 上 以 均匀 角度 的 方式 选取 种 子 点 , 对 高 纬度 区 
域 地 磁场 进行 了 磁力 线 可 视 化 , 并 比较 了 不 同 物理 参 
数 对 地 磁场 形状 的 影响 .根据 真实 地 磁场 的 物理 特 
征 , 在 该 磁 经 圈 上 选取 种 子 点 能 使 得 磁力 线 大 致 分 布 
在 一 个 平面 内 , 有 效 消除 磁力 线 在 视觉 投影 面 的 交叉 
重重 现象 . 但 由 于 是 均匀 角度 方式 选取 种 子 点 , 绘制 
磁力 线 的 下 密 不 能 客观 反映 地 磁场 强度 的 空间 分 布 . 

Hep] 在 磁 赤 道 面 上 距 地 心 距离 为 工 个 地 球 半径 
的 磁 纬 圈 上 以 均匀 角度 的 方式 选取 种 子 点 进行 磁力 
线 绘 制 , 反映 出 了 磁 才 的 基本 轮廓 特征 , 但 由 于 磁 学 
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参数 工 选 取 的 任意 性 , 磁力 线 之 间 没 有 约束 关系 , 同 
样 不 能 从 磁力 线 疏 密 程 度 客 观 反 映 磁 场 强度 的 空间 
分 布 , 并 且 视 觉 空 间 中 的 交叉 重 登 现象 也 比较 严重 . 
2014 年 Zhong 考虑 磁力 线 玻 密 与 磁场 强度 空间 分 
利用 磁 赤 道 面 上 磁力 线 划 密 的 约束 规则 反 推 球面 上 
磁力 线 种 子 点 的 选取 规则 , 对 地 磁 偶 极 子 场 进行 可 视 
化 ， 该 方法 能 够 通过 磁力 线 玖 密 客 观 反 映 地 磁场 强 
度 的 空间 分 布 , 但 文献 [3] 只 基于 偶 极 子 场 模型 进行 
了 验证 , 并 没有 在 更 真实 或 更 复杂 的 地 磁场 模型 下 进 
行 应 用 , 其 适用 性 有 待 进一步 讨论 ，. 
针对 在 磁 经 圈 上 以 均匀 角度 方式 选取 种 子 点 方 
法 的 不 足 , 本 文 设 计 了 磁场 强度 线 积分 等 分 算法 , 利 
用 磁 经 圈 上 磁场 强度 分 布 对 种 子 点 分 布 进行 约束 . 
将 该 算法 应 用 于 对 IGRF 模型 和 T96 模型 地 磁 内 外 
源 场 的 可 视 化 绘制 , 并 对 磁力 线 追 踪 结果 中 出 现 的 宛 
余 磁 力 线 进行 过 滤 , 实现 了 对 该 算法 的 有 效 性 验证 . 


1 磁场 强度 线 积分 等 分 算法 


磁场 强度 线 积分 等 分 算法 的 中 心思 想 是 将 特定 
路 径 (如 人 磁 经 圈 、 磁 纬 圈 ) 上 磁场 强度 的 线 积分 进行 
等 分 . 在 磁场 强度 较 大 区 域 , 相 邻 两 个 等 分 点 之 间 的 
距离 较 短 , 反之 则 较 长 , 因此 利用 磁场 强度 线 积 分 等 
分 点 作为 种 子 点 能 够 有 效 反 映 地 磁场 强度 分 布 . 

为 了 避免 视觉 上 的 交 又 重合 现 象 , 参照 文献 [1]， 
在 与 太阳 风 方 向 平行 的 磁 经 圈 上 选取 种 子 点 . 首先 ， 
求 出 磁 经 圈 上 磁场 强度 的 线 积 分 总 量 , 即 


5 =IBl. dl (1) 


其 中 , B 为 磁场 强度 . 
假设 要 在 该 磁 经 圈 上 选取 N 个 种 子 点 , 则 需要 
求 出 该 磁 经 圈 上 磁场 强度 线 积分 的 NN 等 分 点 的 位 
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置 . 将 每 个 种 子 点 到 地 心 的 连 线 与 地 磁 轴 的 夹 角 记 
为 6 ( 取 值 范 围 为 0 到 27), 记 第 i 个 种 子 点 对 应 
的 夹 角 为 0; (i = 0,1,2,…,N 一 1). 选取 00=0 对 
应 的 点 作为 首 个 种 子 点 , 相 邻 两 个 种 子 点 之 间 的 磁场 
强度 线 积分 均 为 B/N, 从 9 到 4% 的 磁场 强度 线 积 
分 


Oi 
而 =| 1Bldl(b) = 
00 


| IBIRA9G =i, i=0,1,2,. ,Nl1. 
0o N 
(2) 


通常 根据 式 (2) 便 可 求 出 0; 的 值 . 但 在 实际 计 
算 中 , 磁场 强度 数据 不 是 连续 的 , 需要 将 磁 经 圈 按 角 
度 均 匀 方 式 离散 化 为 M 个 采样 点 , 采样 点 到 地 心 的 
连 线 与 地 磁 轴 的 夹 角 


27 
oz = ky k=0,1,2,...,M—1. (3) 


第 0 个 采样 点 沿 磁 经 圈 一 周 并 回 到 第 0 个 采样 
点 对 应 的 周 角 qm = 27, 并 记 ax 对 应 的 采样 点 处 
磁场 强度 为 |Bx|， 简 单 起 见 , 假设 第 及 和 第 十 1 
两 个 相 邻 采样 点 之 间 的 磁场 强度 是 均匀 的 , 用 |Bx| 
和 |Bwti| 的 平均 值 来 表示 . 当 上 不 小 于 1 时 ,第 一 1 
和 处 个 采样 点 之 间 的 磁场 强度 线 积分 可 表达 为 
_ [Brl+|Br-1| 27R 
et 


DB 3 


k=1,2,3,...,M. (4) 


当 =0 时 , 记 84 ==0. 对 于 第 i(i=0,1,2,…,N 一 
1) 个 种 子 点 ， 寻 找 非 负 整数 s, 使 得 0; 处 于 os 
与 as+l 之 间 , 则 s 应 满足 


S S 十 1 
DoSG< YG, s=0,1,2,...,M-1. (5) 
j=0 j=0 


式 (5) 中 的 到 可 以 由 式 (2) 得 到 . 由 于 假设 相 
邻 采样 点 之 间 的 磁场 强度 是 均匀 的 , 找到 第 i 个 种 子 
点 对 应 的 s 后 可 通过 线性 插值 得 到 


sD 
一 -站 
0; 二 Qs 十 (Qs+1 2 Qs), 
s 二 1 
§s=0,1,2,...,M—1. (6) 


将 式 (2)~(4) 代入 式 (6) 即 可 求 出 4 的 值 , 进而 可 求 
出 种 子 点 的 空间 位 置 . 


对 磁场 强度 线 积分 等 分 算法 选取 的 种 子 点 进行 
磁力 线 追 踪 , 会 产生 两 种 磁力 线 , 即 封闭 磁力 线 和 开 
放 磁 力 线 ， 在 封闭 磁力 线 中 可 能 会 产生 磁力 线 距离 
过 近 甚 至 重生 的 现象 , 这 意味 着 两 条 磁力 线 之 间 的 区 
域 磁场 强度 非常 大 , 甚至 趋 于 无 穷 , 与 物理 事实 不 相 
符 , 产生 磁力 线 元 余 . 而 开放 磁力 线 则 不 存在 这 类 问 
题 ， 这 种 现象 是 由 封闭 磁力 线 两 次 穿越 选取 种 子 点 
的 磁 经 圈 而 引起 的 . 根据 磁力 线 上 磁场 方向 的 不 同 ， 
可 将 封闭 磁力 线 区 域 细 分 为 磁力 线 发 出 区 域 和 磁力 
线 射 入 区 域 , 磁力 线 沿 磁场 方向 从 磁力 线 发 出 区 域 出 
发 , 进入 到 磁力 线 射 入 区 域 . 磁力 线 发 出 区 域内 由 种 
子 点 确定 的 磁力 线条 数 与 该 区 域内 磁场 强度 分 布 是 
相 匹 配 的 . 但是, 从 磁力 线 射 入 区 域 的 种 子 点 追踪 出 
来 的 磁力 线 会 进入 到 磁力 线 发 出 区 域 , 使 得 最 终 存在 
于 磁力 线 发 出 区 域 的 磁力 线条 数 大 于 磁力 线 发 出 区 
域 的 种 子 点 数 . 同 理 , 磁力 线 射 入 区 域 也 存在 与 磁力 
线 发 出 区 域 类 似 的 现象. 

理论 上 讲 , 从 磁力 线 发 出 区 域 发 出 的 磁力 线条 数 
与 进入 磁力 线 射 入 区 域 的 磁力 线条 数 是 相等 的 . 对 于 
这 两 组 磁力 线 , 只 需要 一 组 便 能 够 反映 整个 封闭 磁力 
线 区 域 的 磁场 强度 分 布 ， 磁 场 强度 线 积分 等 分 算法 
会 同时 产生 这 两 组 磁力 线 , 于 是 就 产生 了 磁力 线 元 余 
现象 , 需要 对 元 余 的 磁力 线 进 行 过 滤 . 


2 实验 绘制 


地 磁场 可 以 分 为 内 源 场 和 外 源 场 , IGRF 模型 向 
和 T96 模型 是 是 国际 上 广泛 使 用 的 内 外 源 场 模 
型 . 本 文选 用 IGRF 模型 和 T96 作为 地 磁场 数据 源 . 
以 1997 年 1 月 1 日 18:05LT 的 数据 为 例 , 输入 参 
数 为 : 太阳 风 质 子 密度 N -= 4.08cm-3; 太阳 风速 
度 viGsE = 一 414.5km's-1，uwyGSE = —34.5km-.s-!, 
vzGSE 二 一 24.4kms-1; 太阳 风 动 压 Paysn = 1.37nPa; 
行星 际 磁场 强度 Bycsp = 一 2.08nT,，BzGsp = 
0.71nT; 磁 暴 环 电 流 指数 6 SYM-H = -3nT. 

由 太阳 风速 度 三 分 量 可 以 求 得 太阳 风速 度 方向 ， 
在 地 磁 (MAG) 坐标 系 门 中 , 太阳 风速 度 方向 在 磁 
赤道 面 内 的 投影 与 MAG 坐标 系 x 轴 夹 角 即 为 种 子 
点 所 在 平面 的 磁 经 度 . 取 距 地 心 3 个 地 球 半径 的 磁 经 
图 来 进行 种 子 点 选取 , 选取 60 个 种 子 点 . 为 了 呈现 
最 佳 视觉 效果 , 将 种 子 点 空间 坐标 转换 到 GSW 坐标 
系 中 进行 磁力 线 绘制 , 因为 GSW 坐标 系 的 Oxz 平 
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面 与 太阳 风 方 向 基本 平行 , 基于 GSW 坐标 系 的 Ozz 
平面 进行 投影 可 以 有 效 减少 视觉 上 的 交叉 重 闪 现象 . 

确定 种 子 点 空间 位 置 后 便 可 进行 磁力 线 绘制 . 生 
成 磁力 线 的 方法 主要 有 两 种 ， 即 数值 积分 法 四 和 
双流 函数 法 RP- 囊 ， 本 文 使 用 的 是 数值 积分 法 ， 利 
用 GEOPACK (2008) 工具 包 0 中 的 五 阶 龙 格 库 
塔 法 进行 磁力 线 追 踪 . 对 于 同一 个 种 子 点 , 需要 沿 与 
磁场 方向 的 正 反 两 个 方向 分 别 妃 踪 , 如 果 一 条 磁力 线 
追踪 出 来 的 点 到 地 心 距 离 小 于 1 个 地 球 半径 ,或 大 
于 100 个 地 球 半径 , 或 追踪 点 超过 2000 个 点 , 则 停 
止 追 踪 . 

对 于 绘制 结果 中 的 元 余 磁 力 线 ， 需 要 进行 过 滤 
处 理 . 可 以 将 磁力 线 区 分 为 封闭 磁力 线 和 开放 磁力 
线 ， 如 果 追 踪 出 的 磁力 线 两 端点 到 地 心 距离 均 小 于 
等 于 1 个 地 球 半径 ,， 则 标记 为 封闭 磁力 线 ; 如 果 追 
踪 出 的 磁力 线 一 个 端点 的 地 心 距离 小 于 1 个 地 球 半 
径 , 男 一 个 端点 的 地 心 距离 大 于 1 个 地 球 半 径 , 则 标 
记 为 开放 磁力 线 . 保留 所 有 开放 磁力 线 , 过 滤 掉 封闭 
磁力 线 区 域 的 元 余 磁力 线 ， 假 设 要 保留 磁力 线 射 入 
区 域 的 种 子 点 追踪 出 来 的 封闭 磁力 线 , 只 需 判断 该 封 
闭 磁力 线 的 种 子 点 沿 磁场 正方 向 追踪 出 的 第 一 个 磁 
力 线 序列 点 是 远离 地 心 还 是 趋向 地 心 . 如 果 是 远离 地 
心 , 说 明 该 种 子 点 位 于 磁力 线 发 出 区 域 , 则 应 舍 去 该 
种 子 点 对 应 的 磁力 线 ; 如 果 是 趋向 地 心 , 说 明 该 种 子 
点 位 于 磁力 线 射 入 区 域 , 则 应 保留 该 种 子 点 对 应 的 磁 
力 线 . 

使 用 磁场 强度 线 积分 等 分 算法 在 磁 经 圈 上 选 
取 60 个 种 子 点 并 进行 元 余 磁 力 线 过 滤 , 得 到 的 绘制 
效果 如 图 1 所 示 . 在 地 球 附 近 , 磁力 线 基本 呈现 出 偶 
极 磁场 的 结构 特征 , 这 体现 了 地 球 主 磁场 的 特征 . 在 
向 日 侧 , 磁力 线 被 压缩 , 形成 磁 层 顶 ; 在 背 日 侧 , 磁力 
线 被 拉 伸 而 形成 磁 尾 , 并 且 南 北 两 个 极 尖 区 结构 也 很 
清晰 , 体现 了 太阳 风 与 磁 层 相互 作用 的 效果 . 结果 表 
明 , 基于 磁场 强度 线 积分 等 分 算法 得 到 的 可 视 化 结果 
描绘 出 了 地 磁场 的 基本 形态 特征 . 


3 效果 分 析 


3.1 ”磁力 线 疏 密 与 磁场 强度 空间 分 布 的 匹配 性 
为 分 析 绘 制 磁 力 线 朴 密 与 磁场 强度 空间 分 布 的 

匹配 性 , 针对 360 个 种 子 点 的 实例 进行 统计 . 采用 磁 

经 圈 上 每 5° 区 间 (0% 一 5°, 5° 一 10°,……., 355。 360。 
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图 1 选取 60 个 种 子 点 并 过 滤 元 余 磁 力 线 后 的 结果 
Fig.1 Results by selecting 60 seed points and filtering 


redundant magnetic lines 


区 间 ) 内 穿 过 的 磁力 线条 数 描述 磁力 线 芯 密 , 同时 采 
用 该 区 间 中 位 点 (2.5°, 7.5°,.… ,375.5°) 的 磁场 强度 
作为 该 区 间 内 的 平均 磁场 强度 . 

图 2 给 出 了 磁力 线 朴 密 与 地 磁场 强度 空间 分 布 
的 关系 . 从 图 2 可 以 看 出 , 两 极 区 域 磁场 强度 较 大 而 
赤道 区 域 磁场 强度 较 小 , 相应 地 两 磁极 区 域 磁力 线 分 
布 比 较 密 集 而 赤道 区 域 磁力 线 分 布 比 较 稀 玖 , 符合 地 
磁场 的 基本 特征 ， 并 且 磁 力 线 政 密 与 磁场 强度 之 间 
呈现 出 明显 的 正 相关 性 , 这 说 明了 磁场 强度 线 积分 等 
分 算法 的 有 效 性 . 
3.2 ”元 余 磁力 线 过 渡 

为 了 从 视觉 效果 上 清晰 分 辨 元 余 磁 力 线 过 滤 的 
效果 , 此 处 选取 60 个 种 子 点 的 实例 进行 分 析 . 图 3 
和 图 4 均 为 采用 磁场 强度 线 积分 等 分 算法 得 到 的 种 
子 点 进行 磁力 线 追 踪 后 的 结果 . 图 3 给 出 的 是 尚未 
对 元 余 磁力 线 进行 过 滤 的 情况 , 图 4 给 出 的 是 对 元 
余 磁力 线 进行 过 滤 后 的 结果 . 

从 图 3 所 示 的 地 球 周围 结果 可 以 看 出 , 磁力 线 
出 现 了 两 两 成 对 的 现象 , 而 且 越 靠近 地 球 , 这 种 现象 
越 严重 , 其 至 出 现 重 从 , 这 就 是 磁力 线 元 余 现象 . 图 4 
对 元 余 磁 力 线 进行 了 过 滤 , 磁力 线 由 图 3 中 的 60 条 
减少 为 图 4 中 的 41 条 , 磁力 线 元 余 现 象 消失 , 而 随 
着 重生 磁力 线 的 消除 , 磁力 线 对 地 磁场 可 视 化 的 准确 
性 得 到 提高 , 视觉 效果 更 加 清晰 . 


为 使 地 磁场 可 视 化 绘制 的 磁力 线 能 够 有 效 反 映 
地 磁场 物理 特征 ， 提 出 了 一 种 等 分 磁场 强度 线 积 
的 磁力 线 种 子 点 选取 算法 .利用 地 磁场 IGRF 模型 
和 了 96 模型 进行 实验 绘制 , 得 到 的 可 视 化 结果 能 够 


高 广大 等 : 地 磁场 磁力 线 可 视 化 种 子 点 选取 的 磁场 强度 线 积分 等 分 算法 19t 


ml number of magnetic lines 一 一 intensity of geomagnetic field 


number of magnetic lines per 5° interval 
OO 一 ID 上 wm 人 人 了 家 I % 
intensity of geomagnetic field/nT 


图 2 磁力 线 疏 密 与 地 磁场 强度 空间 分 布 的 关系 


Fig.2 Relation between the density of magnetic lines and the distribution of geomagnetic field intensity 


图 3 选取 60 个 种 子 点 并 包含 元 余 磁 力 线 的 结果 


Fig.3 Results derived from 60 seed points including redundant magnetic lines 
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图 4 选取 60 个 种 子 点 并 过 滤 元 余 磁 力 线 的 结果 


Fig.4 Results derived from 60 seed points excepting redundant magnetic lines 


对 磁 层 顶 、 磁 尾 、 磁 际 等 结构 进行 清晰 表达 , 反映 出 算法 可 对 磁力 线 的 空间 分 布 进行 有 效 控制 , 绘制 得 到 
地 磁场 的 基本 特征 . 的 磁力 线 朴 蜜 与 磁场 强度 空间 分 布 有 较 好 的 相关 性 . 
对 磁场 强度 线 积分 等 分 算法 的 统计 分 析 表 明 , 该 ”此 外 ,算法 中 的 元 余 磁 力 线 过 滤 方 法 能 够 有 效 修正 可 
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视 化 结果 中 的 磁力 线 元 余 现象 , 使 地 磁场 磁力 线 可 视 
化 的 物理 准确 性 得 到 提高 , 视觉 效果 更 加 清晰 . 

磁场 强度 线 积 分 等 分 算法 具有 一 定 的 有 效 性 , 但 
还 有 一 些 不 足 . 本 文 在 选取 种 子 点 磁 经 圈 半 径 工时 
具有 一 定 随意 性 , 选择 不 同 的 工 可 能 导致 可 视 化 最 
终结 果 有 所 不 同 , 而 算法 目前 并 没有 考虑 种 子 点 所 在 
磁 经 圈 的 半径 问题 . 

该 算法 的 适应 性 还 可 以 进一步 拓展 , 例如 这 里 仅 
在 一 个 磁 经 圈 上 选取 种 子 点 , 没有 对 地 磁场 整体 空间 
进行 可 视 化 , 仅 表 现 了 地 磁场 的 一 个 侧面 , 信息 表达 
的 完整 性 有 符 深 入 研究 . 但 是 , 着 考 虑 同时 在 多 个 磁 
经 圈 上 选取 种 子 点 , 存在 多 个 磁 经 圈 之 间 互 相约 束 的 
问题 . 如 有 果 同 时 考虑 磁 经 圈 和 磁 纬 圈 , 两 者 之 间 如 何 
互相 约束 , 在 多 个 磁 经 圈 选 取 种 子 点 时 导致 的 视觉 空 
间 中 的 磁力 线 交 叉 重 一 现 象 等 尚 待 深入 分 析 . 


致谢 ”所 使 用 的 GEOPACK-2008 工具 包 和 T96 模型 
函数 由 N A Tsyganenko 教授 提供 . 国家 科技 基础 条 件 平台 
地 球 系统 科学 数据 共享 网 一 空间 科学 数据 共享 平台 提供 了 
数据 支持 . 


参考 文献 


[1] TSYGANENKO N A. A model of the near magnetosphere 
with a dawn-dusk asymmetry 2. Parameterization and fit- 
ting to observations [J].J. Geophys. Res.: Space Phys., 
2002, 107(A8): SMP 10-1-SMP 10-17. DOI:10.1029/2001- 
JA000220 

[2] HE Huan. Study on Key Technologies for Visualization 


of Space Environment [D]. Beijing: University of Chinese 


[11] 


Chin. J. Space Sci.。 空间 科学 学 报 ”2016, 36(2) 
Academy of Sciences, 2009 ( 贺 欢 . 空间 环境 可 视 化 关键 技术 
研究 [D]. 北京 : 中 国 科学 院 研究 生 院 , 2009) 

ZHONG Jia, ZOU Ziming. The algorithm of density vi- 
sualization of magnetic flux lines based on geomagnetic 
dipole field [J]. Prog. Geophys., 2014, 29(6): 2614-2619 
BAI Chunhua, XU Wenyao, KANG Guofa. Earth’s main 
Geophys., 2008，23(4): 


magnetic field model [J]. Prog. 
1045-1057 

XU Wenyao, DU Aimin, BAI Chunhua. Magnetic field 
model of the Earth’s magnetosphere [J]. Prog. Geophys., 
2008, 23(1): 14-24 

WANLISS J A, SHOWALTER K M. High-resolution glob- 
al storm index: Dst versus SYM-H [J]. J. Geophys. Res., 
2006, 111, A02202. DOI:10.1029/2005JA011034 

XU Wenyao. 
space physics and relevant coordinate systems [J]. Prog. 
Geophys., 2007, 21(4): 1043-1060 

LI Qingyang, WANG Nengchao, YI Dayi. Numerical Inte- 
gration [M]. Beijing: Tsinghua University Press, 2001 ( 李 
庆 扬 , 王 能 超 , 易 大 义 .数值 分 析 [M]. 北京 : 清华 大 学 出 版 社 ， 
2001) 

KNIGHT D, MALLINSON G. Visualizing unstructured 
flow data using dual stream functions[J]. IEEE Trans. 
Visual. Comp. Graph., 1996, 2(4): 355-363 
KENWRIGHT D N, MALLINSON G D. A 3-d streamline 


tracking algorithm using dual stream functions [C]//Pro- 


Data organization in geomagnetism and 


ceedings of the 3rd Conference on Visualization’92. 
Boston: IEEE Computer Society Press, 1992: 62-68 
TSYGANENKO N A. Effects of the solar wind conditions 
in the global magnetospheric configurations as deduced 
from data-based field models[C]//International Confer- 
ence on Substorms. Paris: European Space Agency, 1996: 
181 


